Properties of short channel ballistic carbon nanotube transistors with ohmic contacts 
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We present self-consistent, non-equilibrium Green's function calculations of the characteristics 
of short channel carbon nanotube transistors, focusing on the regime of ballistic transport with 
ohmic contacts. We first establish that the band lineup at the contacts is renormalized by charge 
transfer, leading to Schottky contacts for small diameter nanotubes and ohmic contacts for large 
diameter nanotubes, in agreement with recent experiments. For short channel ohmic contact devices, 
source-drain tunneling and drain-induced barrier lowering significantly impact the current-voltage 
characteristics. Furthermore, the ON state conductance shows a temperature dependence, even in 
the absence of phonon scattering or Schottky barriers. This last result also agrees with recently 
reported experimental measurements. 



I. INTRODUCTION 

Since their original fabrication 1,2 , carbon nanotube 
(NT) transistors have seen much improvement in their 
performance 3 , and in the understanding of the science 
that governs their operation. For example, it has recently 
been demonstrated that the metal used to contact the 
NT has a strong influence on the device behavior, with 
reports of ohmic contacts for Pd 4 *^ and Au&, and Schot- 
tky contacts for Ti£. Furthermore, it has been demon- 
strated that scaling of Schottky barrier NT transistors is 
very different from traditional transistors^. These studies 
highlight the danger of using concepts from conventional 
silicon devices and applying them to NT devices. In such 
circumstances, modeling can play a key role in predicting 
device behavior, and in supporting experimental conclu- 
sions. 

While some theoretical and modeling work beyond tra- 
ditional transistor models has been done to study NT 
transistors, these calculations have used various approx- 
imations, either as simplified treatment of the NT elec- 
tronic properties^, or by simplifying the calculation of 
the current at finite bias voltage^. These approxima- 
tions prevent a more detailed and quantitative study of 
the properties of NT devices, particularly under non- 
equilibrium conditions. Here, we present a full self- 
consistent calculation that treats the NT electronic struc- 
ture and the non-equilibrium quantum electron transport 
from an atomistic perspective. This allows the calcula- 
tion of the charge, electrostatic potential and current in 
the device from a quantum transport approach, and in- 
cludes quantum effects such as energy quantization and 
tunneling. It also naturally captures distortions of the 
NT bandstructure due to the spatially-varying electro- 
static potential and is amenable to extensions to include 
scattering with defects, surfaces, phonons, etc. 

The calculations presented here provide an approach 
for evaluating the performance of ultimate NT tran- 
sistors, those with ballistic transport, ohmic contacts 
and nanometer channels (a recent paper presented a 



similar approach 10 , but focused on Schottky barrier 
NT transistors). We show that side contacts between 
NTs and metal are governed by charge transfer and 
band re- alignment, leading to a NT diameter above 
which contacts are ohmic, in agreement with recent 
experiments^ 1 ^. For the ohmic contact devices, short 
channel effects can dominate the device behavior for rel- 
atively large channel to gate length ratios. Furthermore, 
while the temperature dependence of the ON state con- 
ductance in NT devices can be dominated by Schottky 
contacts or phonon scattering, we show that, even in 
the absence of these effects, the conductance can have a 
strong temperature dependence. In fact, we find that the 
ON state conductance decreases with increasing temper- 
ature, in contrast to Schottky barrier NT transistors and 
in agreement with recent experimental measurements on 
Pd-contacted NT transistors 4 . 



II. MODELING APPROACH 

As shown in figures 1 and 2, the system consists of an 
infinitely long single-wall zigzag NT embedded in metal 
at its two ends and coated by a dielectric in the chan- 
nel region. The gate insulator dielectric constant e is 
equal to 3.9, as for Si02, and in this work the dielec- 
tric is wrapped by a cylindrical gate of radius 3 nm in 
the channel region (unless otherwise stated). The nan- 
otube sits in the middle of a smooth cylindrical hole in 
the metal and dielectric, and we take a spacing of s =0.3 
nm between the NT and the metal or dielectric walls. 
The length of each contact A is equal to 4.17 nm, suffi- 
ciently long to converge the electrostatic potential in the 
contacts. 

The electronic properties of the NT are described us- 
ing a tight-binding framework with one n orbital per car- 
bon atom, and a coupling 7 = 2.5 eV— between nearest- 
neighbor atoms. In our calculations, we assume that the 
electronic structure of the NT is not perturbed by atomic 
interactions with the metal or the dielectric. 
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FIG. 1: Sketch of the NT transistor. The NT is embedded 
in metals at its two ends, and in a dielectric in the channel 
region. The dielectric is wrapped by a cylindrical gate of 3 nm 
radius. The separation between the contacts and the central 
dielectric region is to illustrate the structure in the channel; 
in the calculations, the contacts touch the dielectric. 

To calculate the electronic current flowing through 
the NT, we apply the non-equilibrium Green's function 
formalism^* to the NT device, dividing the NT in princi- 
ple layers, with each layer corresponding to a ring of the 
zigzag NT. The device is divided into a scattering region 
(layers 1 to N) connecting the lead regions, which con- 
sist of semi-infinite nanotubes embedded in the source 
and drain contact metals. This approach is illustrated in 
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FIG. 2: Sketch of a cross-section of the nanotube transistor, 
showing the important device dimensions and the principle 
layers used for the quantum transport calculations. 

For this geometry, the current can be obtained from 

1=^ JdEReG< >N+1 (1) 
where G < is the Green's function, indexed according to 



the principle layers in the device, with TV total layers in 
the scattering region. G < is calculated by solving the 
matrix equations 

G< = G fl £<G flt (2) 

and 

G R = [(E-eU)I-H -X* (3) 

where Hq is the tight-binding Hamiltonian for the iso- 
lated NT and U is the electrostatic potential on each layer 
of the system. In our tight-binding representation, the 
Hamiltonian matrix elements are H p,2p ~ 1 — Hq P ~~ 1,2p — 
2 7 cos(^), H 2 p - 2p+1 = H 2p+1 > 2p = 7 where M is the 
number of atoms around a NT ring and J = 1, M la- 
bels each of the NT bands. Because of the cylindrical 
symmetry, we take the electrostatic potential on every 
atom of a ring to be the same, allowing substantial re- 
duction in the size of the matrices (from M N x MN to 
N x N). The functions r and £< represent the inter- 
action of the scattering region with the semi-infinite NT 
left (L) and right (i?) leads. These functions are given 
by 

^LM = 4,r9l, R tl,r (4) 

and 

£< = -2tIm(/iE£ + /B£g) (5) 

where tl^r is the matrix that couples the scattering re- 
gion to the left and right leads, respectively, g E ' R is the 
Green's function for the semi-infinite left and right leads 
and Jl,r are the left and right fermi functions. For the 
zigzag NT in Fig. 2, and within a nearest-neighbor tight- 
binding representation, only the first and last rings cou- 
ple with the leads. Therefore, the coupling matrices T£ jj 
each have only one non-zero element, equal to 7, and we 
have the non-zero components 

[sf] n = i 2 m m (6) 

[ s fl]ivjv = T 2 [#fl]jV+liV+l 

where the indices of g R {Sr) run from —00 to 
(N + 1 to 00). The surface Green's functions g R R are 
calculated using an iterative layer doubling technique^, 
while G R is calculated by using an efficient approach that 
iteratively builds up the Green's function at each layer 
across the device region^. 

The electrostatic potential is calculated by solving 
Poisson's equation with a spatially-dependent dielectric 
constant k 

V-(nVU)=--^-a (7) 

in three-dimensional cylindrical coordinates on a grid 
with the source charge on the NT a, and with boundary 
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conditions at the source, drain, and gate surfaces, and at 
the boundaries of the simulation cell. We model the in- 
terface between the oxide and vacuum with the function 



J 



2«(r) = (s + 1) + (e - 1) tanh [(r - R - s) /£] where £ is 
a length scale on the order of a lattice constant. 
In cylindrical coordinates, the above equation is 
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with the boundary conditions 
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where L is the channel length, Vds is the drain-source 
voltage and V gs is the gate-source voltage. As will be 
discussed below, we use a uniform distribution of the 
charge in the azimuthal direction, thus reducing the 
three-dimensional partial differential equation to a two- 
dimensional partial differential equation in r and z. To 
solve this two-dimensional equation numerically, we use 
a finite-difference scheme 
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where the subscripts i,j denote a function evaluated at 
position (r,, Zj) with rj = (i — l)A r and Zj = (j — 1)A Z , 
and i c — [R + s) / A r is the radial grid point where 
the contacts and dielectric begin. The dielectric/vaccum 
boundary term was evaluated using a forward difference 
scheme, assuming that ( < A r ; in practice we use a value 
£ = 0.0085 nm and we verified that the ratio of the elec- 
tric fields equals the ratio of the dielectric constants at 
the interface. The finite-difference boundary conditions 
are 
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where js = A/ A z is the axial grid point where the source 
contact ends, jr> = {L + A) /A z is the axial grid point 
where the drain contact begins, and ic is the radial grid 
point at the gate radius. In our calculations we use 
grid spacings of A 2 =0.07 nm in the axial direction and 
A r =0.03 nm in the radial direction. Numerical solution 
of Eq. I|10|) is obtained using a successive over-relaxation 



procedure^. Once the electrostatic potential is obtained, 
the value for U on each layer of the system (i.e. that en- 
ters the Hamiltonian) is taken as the value of the electro- 
static potential at the atomic position of each ring along 
the NT. 

To obtain the charge density, we note that the tight- 
binding technique only provides the total charge on each 
layer; in our formalism the total charge on layer I is 
given by J dE Im , which needs to be spatially dis- 
tributed. We assume a uniform distribution of the charge 
in the azimuthal direction, and spatially distribute the 
total charge in the radial and axial directions with a 
Gaussian smearing function. Thus the three-dimensional 
charge density is given by 



a{r,<j),z) = ^2g(z 



zi,r- 



dE Im G% 



(12) 



where g(z — zi,r — R) = (4:n 2 R<7 z <j r ) 
exp [~(z - zi) 2 /2a 2 ] exp [-(r - R) 2 /2a 2 ] with R the 
tube radius, zi the position of ring I, and a z and ay 
the smearing lengths in the axial and radial directions 
respectively (this expression for g is valid when R^$> a r , 
and we use values of a z = 0.14 nm and ay = 0.06 nm). 

Our overall procedure is to start from some initial guess 
for the potential and solve for the Green's functions from 



the system of Eqs ffi-TBfr. calculate the charge distril 
tion in the device from Eq. 1121 and then re-calculate 1 
potential profile across the device based on solving Pc 
son's equation. This process is iterated until both 1 
potential and charge distribution converge; in pract 
this is done using a simple mixing of the charge at ee 
iteration step. 

As discussed in a previous workS, short channel I 
transistors can possess two different operating regimes 
conventional transistor regime, and a Coulomb blocka 
regime at large gate voltages due to channel inversii 
While simplified models exist to describe the Couloi 
blockade regime^, a proper treatment involves includi 
the electron-electron interaction self-consistently in t 
Green's function approach. In this paper, we focus 
the transistor regime and leave the Coulomb blocka 
regime for future work. 

III. PROPERTIES OF CONTACTS 

Before discussing the full current-voltage charactei 
tics of the NT transistor, we begin by discussing the pre 
erties of the contacts for NTs embedded in metals. F 

3 shows a cross-section of the contact, and the asso 
ated bare band lineup. For such a contact, the differer 
between the metal Fermi level Ep and the valence ba 
edge E v before charge transfer is simply given by 

E F - E v = $ WT - $ m + -E g (] , 

where $,„ and $ nt arc the metal and NT workfunctions 
respectively, and E g is the NT band gap. A positive value 
for Ep — E v indicates a Schottky barrier, while a nega- 
tive value indicates an ohmic contact (here we focus on 
the barrier for holes, the relevant quantity for high work- 
function metals used in most experiments; the barrier for 
electrons is E v + E g — Ep and would be relevant for low 
workfunction metals. Because the NT band structure is 
symmetric around the midgap, the results presented be- 
low apply to both cases). Given the relation E g = aj/d 
between bandgap and NT diameter d (a = 0.142 nm is 
the C-C bond length), the band lineup thus depends on 
the NT diameter. The behavior of Eq. Ellis shown in Fig. 

4 for two different values of 4> m — &nt as a function of 
the NT diameter. This simple picture for the band lineup 
is modified due to charge transfer between the metal and 
NT, as we now discuss. 

Because the metal Fermi level is not in line with the 
charge neutrality level in the NT (midgap for an undoped 
NT), charge transfer between the metal and NT occurs, 
as illustrated in Fig. 3a. From an image potential ar- 
gument, this transferred charge creates a charge dipole 
at the nanotube/metal interface, and an associated elec- 
trostatic potential. This electrostatic potential in turn 
shifts the bands, and changes the amount of transferred 
charge. Thus, the charge and potential must be deter- 
mined self-consistently. This behavior can be described 




FIG. 3: Panel (a) shows a sketch of the cross-section of the 
nanotube/metal contact, while panel (b) is the bare band 
lineup. 

using simple models for the charge and the potential (the 
results of the atomistic approach presented in the previ- 
ous section will be discussed further below). The charge 
per atom on the NT can be expressed as 

a = eN J D(E)f(E - E F )dE (14) 

where 

D(E) = a f \ E + eV ^\ (15) 

n R ~t J(E + eV NT ) 2 -(E g /2) 2 

is the density of states, / (E — Ep) is the Fermi function 
and TV = 4/(3v / 3a 2 ) is the atomic areal density. We 
assume a uniform distribution of the charge on the NT. 

For the geometry of Fig. 3a, solution of Poisson's equa- 
tion gives the potential on the NT as 

t/ eR R + s 

eV NT = -a — In — — . (16) 
£o it 

For a given NT, equations 1141 and 1161 can be solved self- 
consistently. The solid lines in Fig. 4 show results of 
such calculations for the room temperature alignment be- 
tween the metal Fermi level and the NT band edge for 
&m — $nt — 1 eV and 0.3 eV. Comparison with the bar- 
rier predicted from Eq. 1131 (dotted lines in the figure) in- 
dicates that the charge transfer significantly changes the 
band line-up. While for the higher workfunction differ- 
ence the contacts are ohmic down to very small diameter 
NTs (Fig. 4a), for the lower workfunction metal (Fig. 
4b), the results indicate a transition between Schottky 
and ohmic behavior at a NT diameter around 1.8 nm. 
This result agrees with recently published experimental 
dataii^. 

The results of this simplified analytical model can 
be compared with calculations from the atomistic, self- 
consistent, Green's function approach introduced in the 
previous section. For this comparison, the electrostatic 
potential on the NT from the atomistic approach is taken 
to be £7i Hi i where in = R/A r is the grid point at the NT 
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FIG. 4: Band lineup at nanotube/metal contacts. Dotted 
line is prediction from Eq. QH^ , solid line is computed from 
Eqs 1141 and 1161 . and circles are calculated from atomistic 
approach. 

radius and we use V s d = Vgs — (since the contact length 
is long enough to screen the electric fields from the drain 
or gate electrodes, the value of Ui Rt % is independent of 
the choice of Vd s and V gs ). The data points in Figure 4 
show results of such calculations for several semiconduct- 
ing zigzag NTs, demonstrating excellent agreement with 
the simplified approach presented above. 



IV. TRANSISTOR CHARACTERISTICS 

Having established the conditions needed to make 
ohmic contacts to NTs, we now discuss the room- 
temperature behavior of NT transistors with such con- 
tacts. For explicit calculations, we take the NT midgap 
as the energy reference level, take $ m — &nt — 1 eV, and 
use a (17,0) zigzag NT. This NT has a radius of 0.66 nm 
and in our tight-binding description, it has a bandgap of 
0.55 eV. 

The current versus gate-source voltage characteristics 
of the NT transistor calculated using the atomistic ap- 
proach of section II are shown in Fig. 5 for two different 



values of the channel length, and for source-drain volt- 
ages of -0.1 and -0.3 V. The general behavior consists of 
an ON state at smaller values of V gs and an OFF state 
at larger values of V gs . These two regimes originate from 
the presence or absence of a barrier blocking the hole 
current, as shown in the band diagrams of Fig. 5b. The 
10 nm channel device shows rather poor characteristics, 
with an ON/OFF ratio of only 1000 and a subthreshold 
swing S = (dlogl/dVgsY 1 = 160 mV/decade. This be- 
havior originates from poor control of the gate over the 
electrostatics in the channel, and from tunneling across 
the hole barrier in the OFF regime. Increasing the chan- 
nel length to 20 nm as in Fig. 5b gives a much better 
S = 69 mV/decade and improved ON/OFF ratio. A 
long-channel transistor with undoped channel is expected 
to be maximally OFF at V gs = Q; it is clear from Fig. 5 
that even the 20 nm device is ON at that gate voltage. 
The inset in Fig. 5b shows the calculated band bending 
for V ga = and Vds = 0, indicating that the channel 
is effectively doped p-type due to a long-ranged charge 
transfer from the contacts^. Thus the NT device is nor- 
mally ON, despite the fact that the channel is undoped. 
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FIG. 5: Current as a function of gate-source voltage for chan- 
nel lengths of (a) 10 nm and (b) 20 nm. The insets show the 
band bending for V gs = (left) and V gs — 1 V (right). 

The output characteristics of the NT transistor are 
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shown in Fig. 6, indicating good saturation of the current 
for the 20 nm channel device, but no saturation for the 
10 nm device. This behavior can be attributed to drain- 
induced barrier lowering (DIBL) as indicated in Fig. 7. 
There, it is shown that for the 10 nm channel device, 
application of a source-drain voltage reduces the barrier 
between the source Fermi level and the middle of the 
channel. A signature of the DIBL effect is also seen in 
Fig. 5a, where there is significant shift between the trans- 
fer characteristics for source-drain voltages of —0.1 and 
—0.3 Volts. A quantitative calculation of the DIBL from 
AI/AVd s at 0.1 \±A gives a very large value of 310 mV/V. 
Fig. 5b shows that the DIBL for the 20 nm channel 
length is reduced substantially to 11.5 mV/V, consistent 
with the saturation of the output characteristics of Fig. 
6. 



we repeated the calculations of the subthreshold swing 
for gate radii of 1.5, 6 and 9 nm; the inset in Fig. 8a 
shows that a good scaling behavior can be obtained if 
the channel length is scaled by a quantity a ~ WR^+l, 
where I = 1 nm. While we have not been able to derive 
an analytical expression to justify this scaling behavior, 
we note that electrostatic analyses for cylindrical gate 
transistorsiSiSS predict a scaling quantity proportional to 
the oxide thickness a ~ (R g — R — s); we have tried this 
type of relation and find that it does collapse the data 
onto a single curve. We suspect that the unusual dielec- 
tric response of the NT—, strong charge transfer from 
the contacts and the actual device geometry render the 
conventional analyses inapplicable; more work is needed 
to address these issues. 
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FIG. 6: Current as a function of drain-source voltage for (a) 
10 nm channel length and (b) 20 nm channel length. 



The reduction of the subthreshold swing and the DIBL 
with increasing channel length is captured in Fig. 8. It 
is clear that both quantities vary rapidly over the same 
length scale of a few nanometers. Typical device con- 
straints require DIBL less than 100 mV/V and S less 
than 80 mV/dec; this restricts the nanotube transistor 
in this geometry and with this dielectric constant to a 
channel length greater than 15 nanometers, outside of 
the shaded areas in the figure. To extend this analysis, 
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FIG. 7: Calculated self-consistent band bending for the 10 nm 
channel device for a gate-source voltage of 0.8 V. The solid 
lines are the valence and conduction band edges; horizontal 
dotted lines are the metal Fermi levels. 
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FIG. 8: Panels (a) and (b) show the variation of DIBL and 
subthreshold swing on channel length, respectively. Shaded 
areas are regions where the short channel effects are larger 
than typical device requirements. The inset in panel (a) shows 
collapse of the data for several devices of different channel 
lengths and gate radii upon scaling of the channel length. 
The horizontal dotted line in panel (b) is the theoretical limit 
for the subthreshold swing. 



the temperature dependence of the conductance C, we 
note that in the ON state the bands are flat with the 
Fermi level in the valence band, and we consider the 
expression 14 



4e 2 

c = - 



(17) 



where P is the transmission probability and / is the 
Fermi function. Assuming that P(E) — 1 for E < E v — 
Ep + A, we obtain 



C = 



4e 2 



s A/kT 



h 1 



„A/kT • 



(18) 



Thus, the conductance is sensitive to the self-consistent 
position of the metal Fermi level with respect to the va- 
lence band edge, A = E v — Ep. In principle, the value of 
A depends on the contact capacitance and temperature, 
and requires numerical evaluation at each temperature. 
But to illustrate the behavior, the dotted lines in Fig. 9 
show the results of Eq. (fTSfl for a constant A equal to 
its zero temperature value (we scaled the expression to 
fit the zero temperature conductance calculated numeri- 
cally). Equation ljl8|l indicates that the conductance can 
vary between the zero temperature limit of Ae 2 /h and the 
high temperature limit of 2e 2 /h, with the variation oc- 
curring over a temperature range kT ~ A. Higher work- 
function metals should thus show a weaker dependence 
on temperature. Indeed, the metal with $ m — $att = 
1 eV gives A as 0.1 eV and its conductance varies more 
slowly than the metal with $ m — &nt — 0.6 eV that 
has A w 0.02 eV. The important point however is that 
the device conductance can vary with temperature by a 
factor of 2 even though there is no phonon scattering or 
Schottky barriers. Figure 9 also shows that recent exper- 
imental data— for NT devices with Pd contacts agree well 
with the predicted behavior for A = 0.02 eV. 



V. TEMPERATURE DEPENDENCE OF ON 
STATE CONDUCTANCE 

As we discussed above, one of the important issues for 
carbon nanotubes is the nature of contacts between NTs 
and metal o 11 ! 12 ! 22 . A useful way to elucidate the proper- 
ties of contacts is to study the temperature dependence 
of the device characteristics^ 2 ^. In NT devices, temper- 
ature is useful to study the role of phonon scattering and 
Schottky barriers, as both of these lead to a tempera- 
ture dependence of the conductance. We now show that 
even for the phonon-free, ballistic-ohmic NT device pre- 
sented here, the conductance can have a strong tempera- 
ture dependence. Figure 9 shows the calculated ON state 
conductance from our non-equilibrium Green's function 
approach as a function of temperature for <& m — $ nt — 
1 eV and <I> m — $nt = 0.6 eV. Clearly, the conductance 
decreases with increasing temperature, with a stronger 
effect for the smaller metal workfunction. To understand 



VI. CONCLUSION 

There are several main conclusions from this work. 
First, for nanotubes embedded in metals, we have shown 
that charge transfer between the metal and the NT can 
renormalize the band lineup, and typically locates the 
Fermi level within a few tens of meV from the valence 
band edge. Second, our non-equilibrium quantum trans- 
port calculations indicate that nano-channel NT transis- 
tors can be dominated by short channel effects due to 
source-drain tunneling and drain-induced barrier lower- 
ing, even at channel lengths where Si transistors show 
good characteristics. Finally, for ohmic contacts, the ON 
state conductance decreases with increasing temperature, 
in contrast to Schottky barrier NT transistors. 
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FIG. 9: Conductance of the nanotube transistor in the ON 
regime as a function of temperature. Black and blue solid lines 
are results of the atomistic calculations for $ m — &nt = 1 
eV and 0.6 eV, respectively. Dotted lines are calculated from 
Eq. (1181 . see text for details. The circles are experimental 
data points extracted from Ref^. 
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